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In this paper, we examine the dynamical behavior of a collapsing star via fully dynamical 
approach. We derive and simplify the complete dynamical equation system for a collapsing 
spherically symmetrical star with perfect fluid, and establish some new constraints on the 
initial and boundary conditions. The mass density near the center can be approximately 
derived in the light-cone coordinate system. The results show that, the singularity can not 
form at the center if the equation of state satisfies some increasing and causal conditions 
and the initial distributions of mass-energy density and velocity are suitable. The dynam- 
ical equations can even automatically remove the weak singularity in the initial data. The 
analysis may give some insights into the structure of a star and the nature of the black hole. 
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I. INTRODUCTION 

In general relativity, we have some typical exact solutions such as the Schwarzschild, Curzon 
and Kerr metrics and some of their extensions to the electrovacuum solutions such as the Reissner- 
Nordstrom and Kerr-Newman metrics [1, 2, 3, 4, 5]. All these solutions have singularity and event 
horizon. Under some common assumptions such as positive energy condition and trapped surface 
for the classical matter field, Penrose and Hawking have proved that the singularity in space-time 
will be certainly generated from some regular states[5]. 

The definition of singularity in curved space-time is a subtle problem. Different from the 
solutions of the ordinary field equations, in general relativity we are free to use any coordinate 
system, the infinity of the metric may be only caused by coordinate system. Besides, the singularity 
of the space-time means the boundary of the manifold, where the normal description for the 
space-time is invalid. A way of detecting singularities within a coordinate system is to find that 
curvature invariants become infinite. These are scalar quantities such as R, R^ V R^ V , R^ ua pR^ VOL ^ , 
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which measure the curvature of space-time and if they become infinite, this is a sure sign that a 
region of space-time cannot be extended [5, 6]. A breakthrough in the understanding of space-time 
singularities was the singularity theorem of Penrose [7], which identified general conditions under 
which a space-time must be geodesic incomplete. This was then generalized to other situations 
by Hawking and others. The singularity theorems are proved by contradiction. Their strength 
is that the hypotheses required are very general, and their weakness is that they give very little 
information about what actually happens dynamically [8]. 

However, as a physical description, the incomplete geodesic is rather difficult for manipulation, 
and it is not a good choice to use the future properties of an evolving manifold. The singular 
curvatures also have a weakness, that is, it can be easily overlooked or confused with coordinate 
singularity. For example, R — > oo is easily overlooked similarly to = 0, (r > 0). We completely 
do not know whether R = or R — > oo on the horizon r = 2m in the Schwarzschild space-time, 
because we can not calculate it from the metric by usual means [9]. The most convenient definition 
may be directly using p — > oo. Although it is equivalent to R — > oo in mathematics, they are 
different in physics. This is because, on the one hand, p — ► oo has manifest physical meaning and 
can not be overlooked. On the other hand, as calculated below, p — > oo occurs before the metric 
itself becomes singular [10], so the space-time can be still treated as a 'place' when the mass-energy 
density approaches infinity [6, 8]. 

For the static spherically symmetrical space-time generated by perfect fluid source 

ds 2 = B(r)dt 2 - A(r)dr 2 - r 2 (d6 2 + sin 2 Qdip 2 ), (1.1) 

we have a number of exact solutionsfll, 12, 13] and generating method[14]. In these solutions 
we usually expressed the mass density p and pressure P as the functions of (A, B) and their 
derivatives. Theoretically, such solutions form the general solution to the problem[12]. However, 
in such expressions the properties of the EOS is quite obscure and concealed, and then most of the 
solutions are probably unrealistic in physics. 

In fact, the realistic static asymptotically flat space-time with spherical symmetry can be gen- 
erally resolved by the following procedure. The dynamics for the space-time with perfect fluid can 
be reduced to the following initial problem of an ordinary differential equation system[10, 15], 

M'{r) = AirGpr 2 , M(0) = 0, (1.2) 
(p + P)(47rGPr 3 + Af) 
P (r) = C ?(r-2Af)r ' ^ = Po ' (0) 

in which P = P(p) is the EOS of the fluid, M{r) is the ADM mass within the ball of radial 



coordinate r, C s = \JP'{p) is the velocity of sound in the fluid. For any given po > we get a 



unique solution. The metric components are given by 

„ / 2MV 1 / f R 2UirGPr 3 + M) \ , . 

A= { 1 --) • B = eXP (-/ r( r -2M) dr )' ( L4) 

where R < oo is the radius of the star. p(r) = P(r) = in the region r > R. 

For any given EOS P = P(p), the solution of (1.2) and (1.3) can be easily solved numerically. 
Similar to the case calculated in [10], we find that if the EOS satisfying the following increasing 
and causal conditions 

1 , [ PoP 1 , (7 > 0), 

0<C S <-, C' s (p)>0, P-J ° H (1.5) 

[ sP, (P ~> oo), 

then all solutions of (1.2) and (1.3) are singularity-free. That is, we always have < p < po and 
R s = 2M{R) < R. The condition 7 > 1 is necessary, which leads to the finite radius of the star 
R < 00 by C s — > 0. In the case 7 = 1, we get R — > 00, and then the space-time is no longer 
asymptotically flat. Or equivalently, finite such gas can not form a static star. This result reveals 
that, if EOS satisfies (1.5) and po has no upper bound, for a star with any given large mass, we 
have regular solutions in equilibrium. The EOS of the matter is decisive for the fate of a star. 

There are some skeptics who proposed special matter models which lead to global solutions 
or the violation of the conditions for singularity theorems[16, 17, 18, 19, 20]. The instability of 
naked singularities in the spherically symmetrical space-time with scalar field gravitating source 
was studied by D. Christodoulouf: ]. There are initial data leading to the formation of naked 
singularities but for generic initial data this does not happen. The structure of the singularity is 
similar to the Schwarzschild case. The above analysis also suggests the formation of singularity in 
space-time should be carefully checked by realistic dynamical models. 

A general dynamical analysis includes the hydrodynamics of the matter source, which is too 
complicated to be achieved. Noticing the fact that the asymmetry of the motion of the star, such as 
the rotation, is always alleviating the singularity, the collapsing of a star with spherical symmetry is 
probably the worst case. So if we can make an effectively discussion on this model, it will be much 
helpful to understand the dynamical behavior of the Einstein equation. This is the motivation of 
the paper. 



II. DYNAMICAL EQUATIONS FOR A COLLAPSING STAR 



The line element in the space-time generated by an evolving star with spherical symmetry is 
generally given by[15] 

ds 2 = u 2 dt 2 - (vdt- wdr) 2 - r 2 {d6 2 + sin 2 9dip 2 ). (2.1) 

Here we take the light velocity c = 1 as unit of speed. For a normal star, (u, v, w) are continuous 
functions of (t, r) with suitable smoothness. The null geodesic along the radius is described by 

(u + v)dt- wdr = 0. (2.2) 

Assume the solution is f(t,r) = C, where C is a constant. Making light-cone coordinate transfor- 
mation 

t = T(f(t,r)), (2.3) 

where T(f) is any smooth function satisfying djTd^f > 0, then we get the line element equivalent 
to (2.1) as follows. 

ds 2 = abdt 2 + 2Vbdtdr - r 2 (d6 2 + sin 2 9dip 2 ), (2.4) 

where (a, b) are continuous functions of (t, r) with suitable smoothness until the star becomes 
singular. (2.4) is similar to the "Eddington - Finkelstein coordinates". The advantage of this 
coordinate system is that, the solution near the stellar center to the Einstein equation can be 
approximately solved under some ordinary conditions. 

For the external Schwarzschild solution, we have conditions 

R 

6=1, a = l -, (for r > R > R s ), (2.5) 

r 

where R = R(t) and R s are respectively the stellar radius and Schwarzschild radius 

f R(t) 

R s = 8irG / p gI ^r 2 dr, (2.6) 
J o 

in which Pg TSuV is the total gravitating mass-energy density including influences of pressure and 
momentum. The detailed definition of p grSuV is given below. 
Denote the 4-vector speed of the fluid by 

U^ = {U, V, 0, 0}, Up = {abU + VbV, VbU, 0, 0}, (2.7) 



which satisfies the line element equation 

1 = g^U^U v = (abU + 2y/bV^ U. (2.8) 

Solve U from it, we get 



Vb ( Va + V 2 + v) 



l - T {^/a + V 2 -V). (2.9) 



U 



So between (U,V), only V is independent variable. When \V\ <C 1, V is approximately equal to 
the usual speed. 

For the perfect fluid model, the nonzero components of the energy-momentum tensor = 
(p + P)U^U V - Pg^ u are given by 

T tt = b(p + P) (aVbU + v) 2 - abP, (2.10) 

Ttr 



b(p + P)\ay/bU + V)U -VbP = T rt , (2.11) 
T rr = b(p + P) U 2 , T ee = Pr 2 , T w = Pr 2 sin 2 6, (2.12) 

in which P = P{p) is the given EOS. The reasonable EOS should satisfy the following conditions 

0<C S = ^<^, (p>0), (2.13) 

where C s is the velocity of sound in the fluid. For the following discussions, we only need the 
condition < C s < Co < 1. 

The nonzero components of Einstein tensor are given by 



-(Vbd t 



G tt = — [Vbdta - abd r a) - -^ab(l - a), (2.14) 

G tr = -Vbd r a - -zVbCL - a) = G rt , G rr = — -, (2.15) 
r r z ra 

_ fa ( d r b d r a\ 1 — a 1Z\ 2 _ . 2 n , n 

Gee = \r\2b + ~)~^ 2 ~ + 1 l) ' W = 66 ' ( } 

where the scalar curvature 1Z depends on the second order derivatives of the metric functions (6, a). 

But it is not used in the following discussion, because the related equations are not independent 

of the used ones. 

By detailed calculations, we find only the following 3 equations are independent ones in the 
Einstein equation G^ u = —8irGT^ v , 

d r b = 8TiGr{p + P)b 2 U 2 , (2.17) 



d t a = 8nG(p + P)rV^b(a + V 2 ), (2.18) 
8 r a = -4irGr((p-P) + {p + P)abU 2 ) + ^—^. (2.19) 



By (2.5), (2.6) and (2.18), we learn R s is conserved in an evolutional star. 

Among the energy-momentum conservation law = 0, calculation shows that, only the 
equation U^T^ = is independent, so we get 

U^p + (p + P)U^ = 0. (2.20) 

Simplifying (2.20) and the consistent equation of (2.18) and (2.19) dt r a = d r ta, we get the dynamical 
equation for (p, V), which is a first order hyperbolic differential equation system, 

(1 - C 2 )VbU 2 d tP + \ [2 - (1 - C 2 S ){1 + abU 2 )] d rP + 2 ^ + + P a l^ U 9rV 

(l-a + 8irGr 2 P) + 2VbUV) (p + P), (2.21) 



1 / all 2 



1 + abU 2 



(1 - C 2 )VbU 2 d t V + \[2 - (1 - C 2 )(l + abU 2 )] d r V + ^p^ rP 

= -— (VbU(l - a + 8irGr 2 P) + 2C S 2 (1 + abU 2 )V^j . (2.22) 
The characteristic speeds are given by 



1 " u + (i - c s )Vbu 2 ~ u(i-c s ) ' {2 - 26) 

v - v - - Cs = v-Cs^Tv 2 - 

2 u (i + c s )Vbu 2 UQ. + C.) • { ■ ' 

The disturbance of the solution (p, V) propagates at such speed % = Vk- 

Among the equations (2.17)-(2.19) and (2.21)-(2.22), only 4 equations are independent, and the 
redundant equation can be treated as identity. For a normal star, the variables have the following 
range of value, 

0<6<1, 0<a<l, 0<p<oo, |V|<oo. (2.25) 

By (2.17) and (2.19), we learn b(-,r) G C°([0, R]) is a monotonic increasing function of r, but 
a(-,r) G C 1 ([0, oo)) has a positive minimum a m i n > 0. The reasonable boundary condition in the 
asymptotic flat space-time is given by 

at r = : a = 1, d r a = 0, V = 0; (2.26) 

R 

at r = R: 6 = 1, a = l--^, p = 0. (2.27) 

R 

The initial data of (p,V) should satisfy (2.17), (2.19) and the above boundary condition (2.27). 
Besides the initial data should still satisfy a constraint similar to 

SttG f r 

p(0,r)r 2 dr < 1, (2.28) 



r 



n 



or in ordinary function 

8vrGp(0, r)r 2 < 1. (2.29) 

(2.28) or (2.29) means that we can not preset a black hole or singularity inside a star, because we 
can not use singular conditions to analyze singularity. The exact form of (2.29) will be derived 
below. 

The Einstein equation has a scaling invariance[22]. Making a transformation similar to that in 
[10], we can get the dimensionless form of equation system (2.17)-(2.19) and (2.21)-(2.22). 
Define an auxiliary function by 

F = ( p + P )abU 2 = {p + P) V ^1Y1~1 . (2.30) 

For a static star, we have F = p + P. Substituting (2.30) into (2.19), we get the solution as 

a = i _ f(p _ P + F)r 2 dr. (2.31) 

r Jo 

In contrast (2.31) with (2.5), we get the total gravitational mass-energy density defined in (2.6), 



(2.32) shows the difference between p gra v and the usual mass density p. By (2.26) and (2.27), we 
find that F satisfies 

F(t,0) = p(t,Q) +P(t,Q), F(t,R) = 0. (2.33) 

III. SIMPLIFICATION OF THE EQUATIONS 

The above equations have a weakness, that is, the geometrical variables (a, b) and mechanical 
variables (p, V) couple each other in a complicated manner, which increases the difficulties for 
discussion. For example, we can hardly treat F as a mechanical variable or a geometrical one. 

The problem can be solved by making the following transformation, 



U = -=, V = V^smhv, (3.1) 
V ab 



where v(t, 0) = 0, then we have 



1/^= (Vab cosh v,-=, 0,0] . (3.2) 



The energy function (2.30) reduces to the pure mechanical quantity 

F = (p + P)e- 2v . (3.3) 

The dynamical eqaution system (2.17)-(2.19) and (2.21)-(2.22) becomes 

r b = 8TrGr{p + P)a- 1 be~ 2v , (3.4) 

8 t a = i7rGraVb(p + P)sinh(2v), (3.5) 

d r a = 8nGre~ v (P sinhv — pcoshv) -\ , (3-6) 

r 

{1-C 2 )e~ v ( Cl cosh-y + sinh-i; 1 . , 

-at/9 + a I a r p + e o r v (coshw — 3 sinh 



Vb(p + P) ^ ' V P + ^ " 2r 

= 47rGr(psinh-y — Pcosh-y) — — (coshu + sinhu), (3.7) 

(l-C s 2 )e-\ a(e v C 2 s d r p ,~ , . , NO 1. - . , , , 

-atv + - I — I- (G s cosh d + smhv)o r v H — (3G S smhw — cosht;, 



v/& 2 V P + P r 



4irGr(C 2 p sinh v — P cosh u) (G? sinh v + cosh u). (3. 

2r 



The characteristic speeds become 



Vi = (C s cosh v + sinh «), (3.9) 

1 — G s 

y 2 = (sinh u - C s cosh t;) . (3.10) 

1 + Gs 

From the above equations, we find (a, b) are separated with (p,P,v). This feature is much 
helpful for theoretical analysis and numeric calculation. Solving (3.4) and (3.6), we get 

a = 1 - - f 87rG> grav r 2 dr, (3.11) 



o 



6 = exp^-8vrG^ (p + P)a- X e~ 2v rdr\ . (3.12) 

In which the total gravitating mass-energy density Pg rscv and the gravitating mass (ADM mass) 
Madm of the star are given by 

Pgrav = e~ v (pcoshv — Psinht;), (3.13) 
f R 

Madm = / 47re v (p coshv - P sinh v)r 2 dr. (3.14) 
Jo 

They are manifestly dependent of (P,v). Madm is conserved due to (3.5). 

Theoretically, we can substitute (3.11) and (3.12) into (3.7) and (3.8), then we get a closed 
equation system for (p,v), which includes all information of an evolving star with spherical sym- 
metry. 



Furthermore, making transformation 



and substituting it into the dynamical equations, we get the coefficients of the equations in the 
polynomial form of w. The physical meaning of w is approximately the usual velocity 



In what follows, we show that, under suitable initial conditions the singularity of the space-time 
can not take place. When (b > 0,a > 0), The space-time is normal and measurable. By (3.5), we 
get 



In which a(0,r) is determined by the initial distribution p(0, r) grav via (3.11). So b(t, 0) — > +0 
certainly occurs before a m i n — > +0. According to (3.12) and v(t,0) = 0, a(t, 0) = 1, we learn that, 
if b(t, 0) — > then we certainly have p(t, 0) — ► oo. That is to say, p(t, 0) — ► oo is a necessary 
condition for the singularity of space-time. 

Since a > 0, by (3.12) we learn < b < 1 is a monotonic increasing function of r. By (2.13), 
we have P oc p as p — > oo. So according to (3.12), b(t,0) — > +0 only if p(t,r) > pQ(t)r~ 2 , (r — ► 0). 
Even if p(t,r) — > po(t)r~ n , (0 < n < 2), we still have 6 > 0, and the space-time is still measurable, 
although the curvature is singular. The following calculation shows that, the dynamical equations 
(2.21) and (2.22) even can automatically remove such weak singularity. So to discuss the singularity 
of the space-time is transformed into to check the order of p(t, r) — > oo as r — ► 0. That is to say, 
we only need to examine the asymptotic behavior of (p, V) near the center of the star while 
p — ► oo(r — > 0). 

The disturbance of (p,V) has only finite propagating velocity (2.23) and (2.24). Within any 
given finite time < t < Tq, the mass density p(t, 0) only depends on the initial distribution of 
(p, V) near the center. The worse case is that, the star has already had extremely high temperature 
and pressure near the center, but it is still collapsing. Here V < means collapse, namely the fluid 
flow towards the center of the star. 

In the state with extreme high temperature and pressure, the EOS of the fluid becomes simple. 
We have the following approximations, 



IV. ANALYSIS OF THE SINGULARITY 




(4.1) 



p=cl 



C S =C < 1. 



(4.2) 



Noticing V(t,0) = 0, by (2.9), (3.11) and (3.12), omitting 0(r 2 ) terms we get 

\V\ « 1, U=^X, (4.3) 

b=b (t), a=l. (4.4) 

Substituting (4.2)-(4.4) into (2.21) and (2.22), we have the simplified dynamical equations 

ITU T + S%T + {B ' V + V> + 4 * GC »> = °' < 4 ' 5 > 

i r^ 8 - l/ + TtW— + (*V + — ) + 4^ = 0, (4.6) 
O 1 + G p r 

where dr = \/b~Qdt. Keep in mind that we have not made approximation for p, so it is valid for all 

range value of p. 

By (4.5) and (4.6), we get a linear first order differential equation for p -1 , 

drp- 1 - drp- 1 + (1 + C^ 2 )d T Vp- 1 + 4vrG(l + Cl)r = 0. (4.7) 

The solution of (4.7) can be exactly solved. This is the advantage of the coordinate system (2.4). 
The solution is given by 

P(r, r) = Pi ° f ^ r) r) exp ((1 + C ~ 2 ) [V(r, r) - V(0, r + r)]) , (4.8) 

in which 

f(r, r ) = l-47rG(l + C 2 )p(0,r + r)- 

/ (T + r-x)exp((l + CQ 2 )[V(x,T + r-x)-V(0,T + r)])dx. (4.9) 
J o 

Since V is a variable of little absolute value, so p(r, 0) — ► oo if and only if /(r, 0) — > 0. 
For a collapsing star, we have V(x,r) < 0, then we have the following estimation, 

exp((l + C - 2 )[y(x,r-x)-y(0,r)]) <exp(-(l + C - 2 )y(0,r)). (4.10) 

Substituting it into (4.9) we get 

f(r, 0) > 1 - 2nG(l + C >(0, r)r 2 exp (-(1 + C^ 2 )V(0, r)) . (4.11) 

Therefore, only if the initial distribution satisfies 

2vrG(l + C 2 )p(0, r)r 2 exp [-(1 + C ~ 2 ) V(0, r)] < 1, p(0, 0) < p < oo, (4.12) 

where po is a constant, the singularity p(r, 0) — ► oo will not occur. (4.12) is the exact form of the 
condition (2.29) mentioned above. 



By (4.8), we find even the initial distribution p(0, r) — ► por~ n , (0 < n < 2), it will become 
regular sonn at r > 0. This is a negative feedback process, so b(t, 0) — > is impossible if the initial 
data is suitable. 

For the ultra-relativistic gas in equilibrium 

V = 0, 8 t p = 0, P=\p, (4-13) 

solving (2.17)-(2.22) we get one solution with singular mass density [15], 

3 4r + 3 

b = a(4r + f3), a = ^-^-, (4.14) 



56irGr 2 ' 7r 
where (a, 3) are constant. Apparently, this solution is nonphysical. If take p(t, 0) = po < °° 5 by 
(1.2) and (1.3), we can check that R — > oo and the solution is not asymptotically flat. So it is 
necessary for 7 > 1 in (1.5) to get a finite star. 
For general critical initial distribution 

where po > is a constant. When r — ► 0, we have (a — > 1, i> — > 0). For any given tq satisfying 
< r - r < ii, by (3.12), we get 

6(0, r) = exp^-87rG(^ + j'J ){p + P) a - l e- 2v rdr^j 

- Ai(ro)exp(-87rG^ r0 ^r ( ir^ 

= ^r^T" 11 ^0, (r -» 0), (4.16) 

in which all < ^4^ < 00 are constants. (4.16) means the space-time itself becomes singular, so 
the condition (4.12) can not be modified. 



V. DISCUSSION AND CONCLUSION 



The above calculation shows that, the property of the EOS of matter source is decisive to the 
fate of the evolving star. If the EOS satisfies the increasing and causal conditions (1.5), which is 
quite normal for ordinary particles, the singularity might be avoided. 

The equation of motion for the matter source in a star is automatically derived from the energy- 
momentum conservation law, so we can not arbitrarily set the mass-energy distribution. Whether 
sufficient matter can be collected in a little ball to form a trapped surface is the key for the reality 
of the singularity theorems. 



As analyzed in [9, 10] and shown above, we find the center of a star is not a balancing point 
for particles, because the gravitational potential converts into the kinetic energy of the particles. 
We can not treat a big star as a solid ball. The concentrated gravitating source might be unable 
to locate at center, it should be distributed at some balancing surface. This phenomenon may be 
the reason why the vacuum solution to the Einstein equation usually has not single singular point, 
which is different from the basic solution of Laplace equation. So for the perfect fluid matter, 
the cosmic censorship principle may be right, although the explanation for the black hole may be 
inadequate [23]. The most external singular surface may be the boundary of valid region of the 
exterior vacuum solution to the metric. 
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